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ABSTRACT 



A digital computer program was developed to find the buckling 
coefficient for rectangular plates with all edges simply supported or 
with all edges clamped. A finite difference technique is used to solve 
the partial differential equation for the deflection of a plate classi- 
cally treated as having only a small deflection compared to its thick- 
ness. The program was prepared to take for an input the stresses at 
nodes formed by grids dividing the plate into rectangles. The stresses 
and deflections at each node are used in writing difference equations. 
An extrapolation formula is featured in the program which allows a 
close approximation to the buckling coefficient without necessitating 
the use of a large number of grid nodes. The program was written in 
FORTRAN 60 but must be used as a FORTRAN 63 program to take advantage 
of some of its inherent flexibilities. Information is provided in 
the output of the program which aids in evaluating the reliability of 



a solution. 
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1 . 



Introduction. 



The determination of the initial buckling load of homogeneous plates 
has been a subject of interest for many years. A very good compilation 
of past approaches may be found in the work by Gerard and Becker [8], 

The methods used were varied and, in general, they may be classified in- 
to two categories: (1) variational methods (commonly known as energy 

methods when applied to mechanics) , the best known of which are those 
attributed to Ritz and B. G. Galerkin and (2) numerical methods in which 
finite difference approximations to the partial differential equations 
of the deflection of the plate at a sufficient number of points on the 
plate result in a set of algebraic equations. The methods under the 
first category result in approximating an infinite set of infinite 
series equations by a finite set of equations which must be solved simul- 
taneously. The second method involves the coefficient matrix of a set 
of algebraic equations whose eigenvalue is found or its determinant 
evaluated for a given value of an appropriate parameter. The parameter 
is varied until the determinant vanishes. 

The first method results in a solution limited to the problem for 
which it was developed. Further, it does not promise to be simple for 
cases which may be more on the practical side. The second method was 
shown to require a considerable amount of labor and time in evaluating 
the determinant of a large matrix, even for a simple case.* Computeriza- 
tion suggests itself as a remedy to such limitations encountered in both 

*A conclusion by J. Yardley in his thesis "Applications of Finite 
Difference Equations to Buckling Problems of Rectangular Plates", 
Washington University, (1948) . 
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